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Non-perturbative renormalization of vector and axial vector currents in 
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Renormalization constants of vector (Zy) and axial-vector (Za) currents are determined non-perturbatively in 
quenched QCD for an RG-improved gauge action and a tadpole-improved clover quark action using the Schrodinger 
functional method. Meson decay constants f p and f n show much better scaling when Zv and Za estimated for 
infinite physical volume are used instead of Z-factors from tadpole-improved one-loop perturbation theory. 



1. Introduction 

In a recent comprehensive study by the CP- 
PACS Collaboration of N f =2 full QCD [T], meson 
decay constants were found to exhibit a very large 
scaling violation over the range of lattice spacing 
a -1 w 1 — 2GeV. This was disappointing since 
an RG-improved gluon action and Sheikoleslami- 
Wohlert quark action with tadpole-improved c sw 
were used. In this calculation, however, one-loop 
perturbative Z-factors, albeit tadpole-improved, 
were used for currents. A natural question 
was whether scaling becomes improved if non- 
perturbative Z-factors are employed instead. 

At Lattice2001, we reported an initial study 
of this problem using the Schrodinger functional 
(SF) method within quenched QCD. We found 
the problem of anomalously large values appear- 
ing in the ensemble of hadron correlators toward 
strong coupling where CP-PACS data of decay 
constants had been taken. In this report, we have 
analyzed this problem in some detail. Here we 
present our final results on the Z-factors includ- 
ing these analysis. 
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2. Method 

We follow the method developed by the AL- 
PHA collaboration [3], and work with a lattice 
geometry of L 3 x T with T = 2L for Zy with a 
vector operator at t = L, and for Za with two 
axial vector operators at t = 3T/8 and t — 5T/8. 

Tree-level values are used for the coefficients of 
boundary counter terms of the action. For im- 
proving the axial current, we adopt the one- loop 
perturbative value for the coefficient ca- 

Values of Zy and Za are determined for f3 — 
2.2 - 8.0 which approximately covers the range of 
the CP-PACS quenched calculation 0, (3 = 2.187 
- 2.575. We have analyzed 200-20000 configura- 
tions depending on (3 value and lattice size. 



3. Exceptional Configurations 

It is straight-forward to calculate Z-factors for 
(3 > 2.6. For lower (3 values on large lattices such 
as 8 3 x 16, however, anomalously large values ap- 
pear in the ensemble of hadron correlators. This 
makes it difficult to determine quark mass pre- 
cisely, and since this means uncertainties in k c , 
also that of Z-t actors. 
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Figure 1. Cutoff dependence of m q at (3 — 2.4 on 
an 8 3 x 16 lattice for three k's around n c . 



We suspect that these "exceptional" configura- 
tions are an artifact of quenched approximation; 
having very small or negative eigenvalues of the 
Wilson-Dirac operator, they would be suppressed 
in full QCD. Since one cannot distinguish "ex- 
ceptional" configurations from "normal" ones on 
some rigorous basis, we restrict the configurations 
used for averaging to those having the value of a 
relevant hadron correlator below some cutoff. We 
then examine if uncertainties under variation of 
the cutoff are contained within some acceptable 
magnitude. 

In Fig. n we illustrate this test for m q for which 
a cutoff is set for fp (see Ref. |3] for definition). 



We estimate k c from 



with the cutoff value 



of 300, because m q is rather stable there. The 
uncertainty in m q at the k c is ~ ±2 x 10 -3 , once 
the cutoff of fp is taken in the range 200 - 1000. 

In Fig. [21 we show how much the Z- factors de- 
pend on m q . Zy is insensitive to m q , and Za 
is consistent within 10% or so, albeit apparently 
exhibiting a more pronounced dependence. 

We analyze the uncertainties in the statistical 
averaging of Z-factors themselves by applying a 
cutoff in f\ , as carried out in Ref. 2 . The conclu- 
sion is similar; Zy is very stable against variation 
of the cutoff, and Za shows a more conspicuous 
variation of 5% or so. 

Uncertainties of Za on an 8 3 lattice of order 
15% in total lead to uncertainty of Za normalized 
at infinite volume of order 30%. The uncertainty, 
however, has little effect in a Pade fit of Za and 
hence final results; Za varies less than 3% at the 
largest coupling (3 = 2.187, even if we artificially 
shift Z A at (3 = 2.4 by 30%. 
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Figure 2. m q dependence of Zy and Z a at (3 
2.4 on an 8 3 x 16 lattice. 
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Figure 3. Size dependence of Zy and Za- 

4. Results for Z- factors 

We determine the Z-factors for infinite volume 
(L* = oo) and also for a fixed finite physical vol- 
ume (L* — 0.8 fm corresponding to 8 3 lattice at 
(3 = 2.6) for comparison. The lattice scale is set 
through the string tension 

As shown in Fig. |3 size dependence of Z- 
factors becomes sizable toward strong couplings. 
Since our quark action employs a tadpole- 
improved value of c sw , we expect 0{a) errors in 
the Z-factors. Therefore we extrapolate or inter- 
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Figure 4. Results for Zy and Z a normalized at 
L* = 0.8 fm and L* = oo. 

polate results linearly in a/L to obtain estimates 
at L* = 0.8 fm and at L* — oo. 

In Fig. 2] we show results of Z-factors as a 
function of bare coupling g 2 , together with Pade 
fits (solid curves in the figure) to them. Non- 
perturbative estimates give values smaller than 
the one-loop perturbative ones (dashed lines) by 
about 20 % (15%) for Z v (Z A ) at the largest cou- 
pling of the CP-PACS simulation, = 2.187. 

Zy PC determined from the ratio of the con- 
served vector current to the local one differs sig- 
nificantly from Zy from the SF method, because 
the local current is not (9(a)-improved. 

5. Scaling Property of Decay Constants 

We compare in Fig. [S] and f p determined 
with non-pcrturbative Z-factors normalized at 
L* = oo (filled circles) with those using perturba- 
tive Z-factors (open up triangles). For compari- 
son, open squares are the results from the stan- 
dard plaquette and Wilson actional]. 

We observe a very encouraging result that with 
the non-perturbative Z-factors scaling violations 
are sizably reduced. Furthermore the continuum 
extrapolation yields values consistent with those 
from the standard action. 

In the same figure, we overlay / w and f p deter- 
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Figure 5. and f p vs. a for our improved action 
with non-perturbative and perturbative (PT) Z- 
factors together with results for the standard ac- 
tion g]. 

mined with Z-factors normalized at finite L* = 

0. 8 fm (open circles). Scaling is best improved 
when Z-f actors are normalized at L* — oo. This 
property is likely related to the fact that 0{a/L) 
errors in Z-factors are removed in the limiting 
procedure L* — > oo. 

We also find that f p determined from the con- 
served vector current (filled down triangle in 
Fig. |5J) exhibits a large scaling violation. 
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